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Abstract. - The asymptotics of the equal-time one-particle Green's function for the half- filled 
one- dimensional Hubbard model is studied at finite temperature. We calculate its correlation 
length by evaluating the largest and the second largest eigenvalues of the Quantum Transfer 
Matrix (QTM). In order to allow for the genuinely fermionic nature of the one-particle Green's 
function, we employ the fermionic formulation of the QTM based on the fermionic ii-operator 
of the Hubbard model. The purely imaginary value of the second largest eigenvalue reflects the 
fcF(= 7r/2) oscillations of the one-particle Green's function at half-filling. By solving numerically 
the Bethe Ansatz equations with Trotter numbers up to TV = 10240, we obtain accurate data 
for the correlation length at finite temperatures down into the very low temperature region. 
The correlation length remains finite even at T = due to the existence of the charge gap. Our 
numerical data confirm Stafford and Millis' conjecture regarding an analytic expression for the 
correlation length at T = 0. 



The one-dimensional (ID) Hubbard model 

1,, 1, 



has been intensively studied in condensed matter physics as one of the fundamental models 
for strongly correlated electrons. In 1968, Lieb and Wu [1] solved the model exactly by means 
of the coordinate Bethe ansatz method. Based on the Bethe ansatz, many physical quantities 
have been studied [2-4]. Especially in 1972, Takahashi introduced the string hypothesis and 
derived the thermodynamic Bethe ansatz equations (TBA) [5], from which several thermal 
properties were calculated for finite temperature [6-9]. 

More recently quite a different approach to the thermodynamics of the ID Hubbard model 
has been developed based on the quantum transfer matrix (QTM) method [10-30]. In this 
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approach, the problem is reduced to the eigenvalue problem of the QTM. The free energy 
and other bulk quantities can be calculated from the largest eigenvalue of the QTM. In 1998, 
Jiittner, Kliimper and Suzuki constructed the QTM utilizing fully the integrability structure 
of the Hubbard model [25] . For the largest eigenvalue they derived a new type of non-linear 
integral equations (NLIE). This set of NLIEs turned out to be very powerful especially for the 
numerical study of the thermal quantities, since it consists of only three auxiliary functions. 

One of the remarkable advantages of the QTM approach is that it enables us to study 
the asymptotics of several correlation functions at finite temperature. Explicit results for 
the correlation lengths can be derived from the ratios of the next-leading eigenvalues to the 
largest eigenvalue of the QTM. Using this idea, Tsunetsugu calculated the correlation length 
of the spin-spin correlations for the half-filled ID Hubbard model by solving the Bethe Ansatz 
equations numerically [18]. In the low temperature limit, his results agree with the finite 
temperature correction of the conformal field theory as oc Vg/T. Later this was confirmed 
analytically by Kliimper and Bariev [22] . 

In this paper, we study another important correlation function of the half-filled ID Hub- 
bard model, namely the correlation length of the equal-time one-particle Green's function 
(c^^Cjo-)- As is well known, the system is an insulator at half-filling and the charge excitation 
has a gap. There is no efi^ect of the charge gap on the spin-spin correlation function, be- 
cause the charge excitation is not relevant for the spin-spin correlations. Quite differently, the 
charge gap plays a significant role for the one-particle Green's function leading to exponential 
asymptotics even at T = 0. The quantitative relation, however, between the charge gap and 
the correlation length is absolutely non-trivial and not well understood yet. In regard to this 
point, Stafford and Millis [31] conjectured that the correlation length ^ for the one-particle 
Green's function at T = is given by 



They obtained formula (2) in a rather heuristic way from an analytic calculation of the finite 
size L dependence of the Drude weight at T = [31] 



Stafford and Millis argued that the correlation length of the one-particle Green's function 
is identical to the length £,c{U) whose computation yielded expression (2). Certainly, in the 
limiting cases [/ ^ and U ^ oo formula (2) shows the expected behaviour [31]. However, 
the status of (2) is that of a conjecture, especially for finite values of U, as a proof is still 
missing for the identity of the finite size scaling length and the Green's function correlation 
length ^. One of the goals of this paper is to confirm formula (2) directly by means of the 
QTM method and numerical calculations. 

In ref. [25], Jiittner et al utilized Shastry's i?-matrix [32-35] to construct the QTM for 
the ID Hubbard model. Therefore, their Bethe ansatz equations are strictly valid only for 
the coupled spin model, which is obtained from the ID Hubbard model through the Jordan- 
Wigner transformation. This transformation poses no problem for the investigation of the 
bulk properties of the system or the asymptotics of correlation functions such as the spin- 
spin correlations. However, for studying genuinely fermionic correlation functions, especially, 
the one-particle Green's function, it is necessary to take properly into account its fermionic 
nature. For this purpose, we employ here the fermionic formulation of the QTM developed 
by one of the authors [36,37]. (A related fermionic QTM for the lattice spinless fermion 
model was studied in ref. [27].) Below we sketch the salient points of the fermionic QTM for 




(2) 



D^{L) = (-l)^/2+iLi/2^(^)e-i/«c(c/)_ 



(3) 
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the ID Hubbard model. First we introduce the fermionic i?-operator for the ID Hubbard 
model [36-38] defined by 



ljk{u,v) = -r!^ {u - v)n<ji^ {u-v) + "° tanh(fe(^) - h{v))n%^ {u + v) 

x7e5H^ + ^^)(2n,-T-l)(2n,i-l), (4) 
T^tki"^) = a{u){-nj„nka- + (1 - - rika)) - b{u){nja{l - rika) + (1 - nj„)nka) 

+c{u){c]^Cka + cl^Cja), {cT =t, i) (5) 



where 

/ \ 7/ \ • / \ 1 smh2h(u) U 
a(u) = cos u, 6(u)=smu, c(u) = 1, = — . (6) 

Due to the Yang-Baxter equation for the fermionic ii-operator [38] 

TZl2 {Ul , U2)^13 {Ul , Uz)TZ2S {U2, U3) = TZ^Z {U2 , 'W3)'^13 (^1 , 'U3)7?.12 {Ul , ^2), 

the row-to-row transfer matrix, t{u) = Str^ {TZaL{u, 0) • • • TZai {u, 0)} constitutes a commuting 

family of operators, [t{u),t{v)] = 0. Hence from the logarithmic derivative of the transfer 
matrix, we obtain the set of commuting local fermionic operators 

t{u) = r(0) exp {ul^') + ^I^^) + +•.•}, (7) 

where I^^^ is nothing but the Hamiltonian H (1). Next, following ref. [36] we introduce 
the super-transposed i?-operator iZajiuj,Ua) = Tlf°{uj,Ua), as well as a conjugated transfer 
matrix f (u) = Stro {1ZaL{0, u) ■ ■ ■ 1Zai{0, u)} . (See the original reference [36] for the definition 
of "st" and the explicit form of TZajiujjUa)-) Since f{u) is actually expressed as 

f(w) = Stra{7^la(0,M)•••7^La(0,^^)}, (8) 

we can easily deduce 

f(u) =T(0)-^exp{-uW + O(u^)}. (9) 
Combining (7) and (9) we obtain the fundamental relation 

Ti-u)fiu) =exp{-2un + 0{u^)} , (10) 

which allows us to represent the partition function of the system as 

Z = Tre-'^^= lim IV {r(-ujv)f (ujv)}^/' 

JV— >oo 



= lim Str{rQ™(«^^,0)}^ Uiv = ^ ) • (H) 

Here /3 is the inverse temperature and N is called the Trotter number. The quantum transfer 
matrix t^'^^{un,v) is defined by 



^QTM 



{un,v) = Trj^^Tla,^j{-UN,v)TZa^_ij{v,UN) ■ ■ ■Tla2j{-UN,v)Tla-,j{v,UN)Y (12) 



From the largest eigenvalue Aq of t^'^^{un, 0), the free energy per site is obtained in the 
thermodynamic limit as 

Pf = - lim InAo. (13) 

JV— >oo 
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The correlation lengths of the (static) correlation functions can be evaluated from the 
next-leading eigenvalues 

i = -5iln(Ai/Ao). (14) 

si 

Though the largest eigenvalue Aq is always real and positive, the next leading eigenvalues are 
not necessarily so. In this case, 91n(Ai/Ao) is responsible for oscillations of the corresponding 
correlation function. 

The fermionic QTM (12) has been diagonalized [36] by means of the quantum inverse 

scattering method based on the technique by Ramos and Martins [39,40]. The eigenstates of 
the QTM are described by the set of rapidities {sj,wi},{j = 1, - ■ ■ ,n,l = 1, - ■ ■ ,m), satisfying 
the Bethe ansatz equations 



X q2{sj-i7) q2{wi+2i-f) qi{wi+ij) 

(p(Sj) = zr^ r^, — -, — ^ = z r^, (15) 

q2{sj+iv g2(w;-2z7) Qiiwi-i-y) 



</.(s) 

with the following settings 

7 = ^, z± = exp(a) {tanujvj'^S (17) 

and 

z{s) = is{l + ^/(l - l/s^)), sinh(a) = - j sin(2uAr). (18) 

Here, for convenience, we have applied the partial particle-hole transformation to the Bethe 
ansatz equations derived in [36]. Note that the one-particle Green's function is "invariant" 

under the partial particlc-holc transformation. 

The corresponding eigenvalue of the fermionic QTM, t'^^^^{un,0) is given by 

n 

A = c{uN)]Jzj, (19) 

where Zj = z{sj). We omit here the precise form of the prefactor c{un), which is in common 
for all eigenstates and is not necessary for the evaluation of the correlation lengths. 

We like to emphasize that the fermionic statistics is properly incorporated in the Bethe 
ansatz equations (15) and the eigenvalue formula (19). For example, we have found that the 
dominant eigenvalue contributing to the correlation length for (c^^Cjo-) is the largest eigenvalue 
in the sector [n = N — l,m = N/2 — 1]. In Fig. 1 (b), we plot an example of the distribution 
pattern of the rapidities Zj for the dominant eigenvalue Aq. As all rapidities Zj turn out to 
be purely imaginary, we plot them in a 90-degree rotated frame, i.e., the horizontal line is the 
imaginary axis. Compared with Zj for the largest eigenvalue Aq (see Fig. 1 (a)), we can roughly 
say that one of the Zj closest to the origin is removed. Note that the complex conjugate of Aq 
is also the largest eigenvalue of the sector [n = — l,m = A^/2 — 1] and the corresponding 
rapidities zj are given by the reflection of those for Aq with respect to the real axis (the vertical 



]J(s-Sj), q2{s) ='[[{s- wi), 



(l-Z-M,s))(l-z+/z(,s)) 
(l + z./z(s))(l + z+/z(s)) 



N/2 



(16) 
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Fig. 1 - (a) The rapidities zj on the imaginary axis (horizontal hnc) for the largest eigenvalue Ao in the 
sector N = 8, n = 8,m = 4, U = 8,(3 = 0.8. (b) The rapidities Zj on the imaginary a^xis (horizontal 
line) for the largest eigenvalue Aq in the sector N = 8, n = 7, m = 3, U = 8,/3 = 0.8 



line in Fig. 1). From the expression (19) and the distribution pattern of Zj in Fig. 1 (b), we 
find Ag to be purely imaginary. This is consistent with the fact that the one-particle Green's 
function should exhibit kp = tt/2 oscillations at half- filling. We have numerically solved the 
Bethe ansatz equations (15) for increasing Trotter numbers N and calculated the correlation 
length ^ from the ratios of the eigenvalues Aq and Aq. The result is plotted in Fig. 2. The 
largest Trotter number we took is N = 10240 down to temperatures as low as T = 0.001. For 
such low temperatures, we have extrapolated our data for = 6144,8192, 10240 to estimate 
the limiting behaviour for oo. The error of the extrapolated results is less than 0.5 %. 

For moderate temperatures (T > 0.01), it was sufficient to use Trotter numbers just about 
N = 4096 to get very accurate data. 

From Fig. 2, we can clearly observe that the correlation lengths remain finite as T ^ 0. 
The extrapolation of our data to T — > is compared with (2) in Fig. 3 (a). To perform the 
extrapolation, we assumed the low temperature behaviour of the correlation length in the form 
^ Ct=o X exp(— 7T) and estimated ^t=o from our numerical data at low temperatures. Ob- 
viously our extrapolations coincide with formula (2) almost perfectly. Therefore we conclude 
Stafford and Millis' conjecture is correct for any value of U. 

Having confirmed that formula (2) gives the correct correlation length ^ at T = 0, we now 
look at its dependence on the gap A, see Fig. 3 (b). Here, we used the well known gap formula 



U J, si 



1 



(20) 



sinh(27rj//?7) 

Usually the correlation length of a gapped system is believed to be simply inversely propor- 



1000 



100 



10 



0.1 



r 








U=1.0 - 




T 


U=2.0 






U=4.0 






U=8.0 - 






Stafford-Millis 


• 











10000 



1000 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
T 




100 



0.02 



Fig. 2 - (a) Temperature dependence of the correlation lengths for one-particle Green's function for 
the half-filled Hubbard model, (b) Similar to (a) at lower temperatures 
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tional to the magnitude of the gap. Equations (2) and (20) indeed give ^ ^ oc A for small U, 
however oc lnA('~ InU) for large U [31]. More precisely, we have 

~A/4 (f/^0), (21) 

^1 ~ln(A/a) (C/^oo), (22) 

where a = [r(i)/-\/27r]^ — 4.37688 (sec ref. [31]). Thus the gap dependence of the correlation 
length docs not follow a simple rule. For a discussion of this issue we like to refer the reader 
to rcf. [41]. 

Finally, for increasing temperature T, the correlation length ^ is getting closer to the 
non-interacting free-fermion {U=0) value, namely ^u=o = 1/ In + ^JJ^TTJTp^'^Pl^ , ir- 

respective of the value of U. Actually, for small C/(~ 1), the correlation length ^ takes almost 
the same value as ^(7=0 even at rather low temperatures (T ~ 0.01) 

In this paper we have presented a numerical treatment of the one-particle correlation length 
^(T) for finite Trotter number N . In a future publication, we will derive the NLIE for the 
corresponding eigenvalue Aq in the limit N ^ 00. This will allow us to obtain analytic results 
for the T ^ limit and to prove formula (2) analytically. 

After completing this manuscript we found a derivation of the zero temperature correlation 
length (2) on the basis of the standard row-to-row transfer matrix of the Hubbard model and 
the one-particle energy- momentum excitation following the method described in [41]. Details 
of this calculation will be published elsewhere. 




Fig. 3 - (a) Comparison of data obtained from an extrapolation of the finite temperature correlation 
length to T = with Stafford and Millis' conjecture, (b) The dependence of the T — correlation 
length ^ on the charge gap A. The dotted hnes correspond to the asymptotics (21) and (22). 
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